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ABSTRACT 

We present results for J and K s near-IR imaging data on a large sample of 88 galaxies 
drawn from the catalogue of Impey et al. (1996). The galaxies span a wide range in 
optical and IR surface brightness and morphology (although they were drawn from a 
catalog constructed to identify low surface brightness galaxies, hereafter LSBGs). They 
were also selected to include very low and high HI mass galaxies in order to ensure that 
they span a wide range of evolutionary states. 
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The near-IR data unveils many features of LSBGs not seen before in the optical. 
First, a high fraction of the observed LSBGs are very luminous in the near-IR, indicating 
that they have a well developed old stellar population, and that older LSBGs are more 
frequent in the universe than data from optical bands suggested. Second, the near-IR 
morphologies are often quite different than seen in the optical. Many diffuse LSBGs 
that are apparently bulgeless when observed in blue bands, instead exhibit nuclei in 
J and K s bands. Third, we find significant trends between the near-IR morphologies 
of the galaxies and their ratio of HI mass to near-IR luminosity. Fourth, we find no 
trend in disk surface brightness with absolute magnitude, but significant correlations 
when the bulge surface brightness is used. Finally, we find that the formation of a bulge 
requires a galaxy to have a total baryonic mass above ~1O 1O M . 

A wide variety of other correlations are explored for the sample. We consider corre- 
lations among morphologies, surface brightnesses, near-IR colors, absolute magnitudes, 
and HI masses. In addition, using previous results by Bell and de Jong (2001), we con- 
vert the galaxies' near-IR luminosities to stellar masses, using color-dependent stellar 
mass-to-light ratios. This allows us to consider correlations among more fundamental 
physical quantities such as the HI mass, the stellar mass, the total baryonic mass, the 
gas mass fraction, the mass surface density, and the metallicity (via the highly metal 
sensitive color index J — K s ). 

We find that the strongest of our correlations are with the ratio of HI mass to total 
baryonic mass, Mhi /Mbaryonio which tracks the evolutionary state of the galaxies as 
they convert gas into stars, and which ranges from 0.05 up to nearly 1 for the galaxies 
in our sample. We find strong systematic trends in how the metallicity sensitive J — K s 
color becomes redder with decreasing Mhi /Mb aryonic , as would be expected for "closed 
box" models of chemical enrichment. However, the increased scatter with increasing 
gas mass fraction and decreasing galaxy mass suggests that gas infall is increasingly 
significant in the gas rich lower mass systems. We argue that the overall range in 
J — K s color argues for at least a factor of 20 change in the mean stellar metallicity 
across the mass range spanned by our sample. We also see strong trends between 
Mni/Mbaryonic and central surface density, suggesting that increased star formation 
efficiency with increasing gas surface density strongly drives the conversion of gas into 
stars. 

Subject headings: galaxies: bulges — galaxies: evolution — galaxies: fundamental 
parameters — galaxies: irregular — galaxies: peculiar — galaxies: photometry - 
galaxies: stellar content — galaxies: structure — infrared: galaxies 
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1. Introduction 

Low surface brightness galaxies (LSBGs) have been a subject of increasing interest in the 
last two decades. Although an initial study by Freeman (1970) suggested that the central surface 
brightnesses of disk galaxies in the Hubble sequence (Sa-Sb-Sc) fell within a rather narrow range, 
Disney (1976) pointed out that the lack of low surface brightness galaxies may in fact be a selection 
effect due to the difficulty in discovering galaxies of very low surface brightness, as had previously 
been recognized by Zwicky (1957). Indeed, many surveys have since turned up large numbers 
of LSB galaxies. In practice, "low surface brightness" has come to mean galaxies whose central 
surface brightness is fainter than 22.0 mag arcsec -2 in the B band-pass (i.e. more than la outside 
of the narrow range defined by Freeman (1970) of /iq = 21.7 ± 0.3-B-mag/arcsec 2 ). However, while 
many galaxies have been found below the "Freeman value" of surface brightness, disk galaxies with 
surface brightnesses that are significantly higher than those of normal spirals do not appear to 
exist. This is shown most clearly in Courteau (1996), who find that there is a rather well-defined 
upper cut-off at a central surface brightness in R of ~ 20.08 mag arcsec -2 . 

It was originally thought that all galaxies with a low surface brightness were early or late- 
type dwarfs (van den Bergh 1959). However, radial velocity observations by Fisher and Tully 
(1975) showed that some galaxies with low surface brightness are actually quite large and luminous 
(McGaugh 1992; Bothun, Impey and McGaugh 1997), for example LSBGs like Malin 1 type. 
Indeed, most surveys have revealed LSBGs with disk scale lengths comparable to the Milky Way 
(although these angular diameter limited surveys are strongly biased toward finding the physically 
largest galaxies). 

Based upon surveys to date, LSBGs do seem to have different properties than those of brighter 
galaxies (possibly indicating different evolutionary histories). The observed trends suggest that 
LSBGs are in large part a continuous extension of the population of normal disk galaxies, reaching 
to higher angular momenta and lower masses (Dalcanton, Spergel and Summers 1997). Previ- 
ous studies (Longmore et al. 1982) showed that LSBGs are much more gas rich (in terms of 
Mg as /M s tars) and bluer than "normal" late type galaxies (McGaugh 1992; de Blok, van der Hulst 
and Bothun 1995; Sprayberry et al. 1995). When combined with their observed low gaseous 
metallicity (McGaugh 1994; de Blok & McGaugh 1997), the evidence suggest that LSBGs have 
had a historically low star formation rate and are relatively unevolved systems compared to the 
bulk of normal spirals. HI observations (van der Hulst et al. 1993; de Blok, McGaugh and van der 
Hulst 1996) show that LSBGs have extended disks with low gas surface densities and high M/L 
ratios, further confirming the unevolved nature of LSBGs. Even external influences seem to have 
had little effect in speeding up the evolution of LSBGs (McGaugh and de Blok 1997). 

Despite the tremendous progress in this field, there are still significant uncertainties in many 
of the above results. First, while the data suggests a young stellar age for LSBGs overall, the 
spread in the integrated optical colors of LSBGs is large, suggesting that they may actually have 
traveled along diverse evolutionary paths. In particular, because of the low surface brightness of 
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the underlying population, only a small fraction of the total number of stars is needed to make the 
colors significantly blue, and thus a small change in the star formation rate can drastically change 
the apparent mean age of the stellar population as inferred from optical observations. Second, 
while emission line measurements have revealed the metallicity of the current gas phase, they have 
not measured the metallicity of the underlying stellar population. Given the large gas reservoirs 
observed in LSBGs and the possibility of substantial late time gas infall, it is not clear that the 
current gas phase metallicity is closely related to the actual metallicity of the stars. Third, the total 
stellar mass in these galaxies (and their contribution to dynamics) is rather uncertain, given the 
large uncertainties in the stellar mass-to-light ratio at optical wavelengths. Finally, recent modeling 
of LSBGs predict that from the assumption that blue LSBGs are currently undergoing a period 
of enhanced star formation, there should exist a population of red, non-bursting, quiescent LSBGs 
(Gerritsen and de Blok 1998; Beijersbergen, de Blok and van der Hulst 1999; Bell et al. 2000; 
van den Hoek et al. 2000). These galaxies should then also be metal-poor and gas-rich, and share 
many of the properties of the LSBGs observed in the optical. 

In all of these areas, we can greatly improve the observational constraints by undertaking a 
systematic study of the properties of LSBGs and normal spirals in the near-IR. Combined near-IR 
and (eventually) optical observations are a perfect probe to help to clarify the true stellar masses 
of LSBGs, to disentangle the roles played by the age and metallicity, and to understand their place 
in the continuum with normal high surface brightness (HSB) spirals. 

We expect that contrary to initial impressions, LSBGs should be readily detectable in the near- 
IR. Although much of the data suggests that the mean stellar population of LSBGs are probably 
young, we know from the very red integrated colors of some LSBGs (such as the giant Malin-type 
objects) that many LSBGs contain a significant old population of stars, and were therefore not 
formed recently. This old population should dominate the near-IR (which is a good tracer of red 
giant and lower mass stars) and thus should be detectable with NIR observations. We also should 
be able to easily identify LSBGs with bulges, as these are among the reddest members of the LSBG 
population McGaugh (1992). There are also indications that the most massive LSBGs, such as the 
super-giant galaxies Malin 1 and Malin 2, are redder than the majority of LSBGs, with luminosities 
that are comparable to those of normal spirals. Thus, near-IR observations should be particularly 
sensitive to the most massive end of the LSBG population. 

While valuable, previous studies of the IR properties of galaxies in the low surface brightness 
regime have been limited due to the difficulty of the observations. Knezek and Wroten (1994) 
observed a preliminary sample of LSBGs in J, H and K, but the sample was very small and biased 
toward massive galaxies. Bergvall et al. (1999) have observed 14 blue LSBGs in J, H and K and 
combined with optical photometry, arriving to the result that many properties of LSBGs observed 
in the optical are reproduced in the near-IR. In particular, they observe the same morphologies in 
B and J, measuring weak optical/near- IR color gradients, a fact which they interpret as a low dust 
content. 
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Near-IR observations have also been used to probe the stellar population differences between 
LSBGs and normal spirals. Bell et al. (2000) undertook a combined optical and IR study of 
26 LSBGs, including 5 from his earlier thesis work in Bell et al. (1999), and examples of both 
normal blue LSBGs and the rarer class of red LSBGs. They compare LSBG average ages and 
metallicities with their physical parameters and find that blue LSBGs are well described by models 
with low, roughly constant star formation rates, whereas red LSBGs are better described by a "faded 
disc" scenario. Among the larger population, strong correlations are seen between an LSBG's star 
formation history and its i^T-band surface brightness, if-band absolute magnitude and gas fraction. 
These correlations are consistent with a scenario in which the star formation history of an LSBG 
primarily correlates with its surface density and its metallicity correlates with both its mass and 
its surface density. Another larger study involving IR data of spiral galaxies was made by de Jong 
(1996), but included only a very small number of true LSBGs. 

The lack of much previous work, and the availability of relatively large IR array detectors, 
both open the possibility to systematically observe a large sample of LSBGs and normal spirals in 
the near-IR, to derive fundamental parameters related to the origin of these galaxies and to their 
relation to normal galaxies. Based upon the large expected rewards and the limited existing data, 
a systematic near-IR study of LSBGs is needed in order to complement existing optical data. We 
undertake such a study in this paper, which is organized as follows. In §2 we present the sample 
used in this work. §3 briefly describes our observations, data reductions, photometric calibrations, 
galaxy photometry, and surface photometry. In §4 we analyze and discuss the most important 
results, including galaxy photometry, central surface brightnesses and scale lengths, the relationship 
between Hi-mass and near-IR luminosities, and the relationship between stellar mass-luminosity 
ratios and baryonic mass. Finally, we conclude in §5. 



2. Sample 

All galaxies were initially selected from the catalogue of large angular size low surface brightness 
galaxies compiled by Impey et al. (1996) from the APM plates, consisting of 516 LSBGs with 
effective radius larger than ~ 10 arcsec. Some constraints were applied to optimize the galaxies' 
observability and their intrinsic stellar and gas content. First, a surface brightness cut was applied 
to ensure detectability at the 40-inch Swope telescope and at the 100-inch du Pont telescope, both 
located at Las Campanas Observatory (LCO). The corresponding central surface brightness cut-off 
was set to Ho{B) = 22.0 mag arcsec^ 2 for the galaxies observed with the 40-inch telescope, and 
to Hq{B) = 23.5 mag arcsec -2 for the galaxies observed with the 100-inch telescope. These limits 
are set considering (1) the corresponding limiting K s magnitudes for the 40-inch and the 100-inch 
telescopes both equipped with a NICMOS3 HgCdTe array: 19.0 mag and 20.5 mag, respectively, and 
(2) the observed average color < B — K >~ 3.0 for LSBGs (Tully and Verheijen 1997). We applied 
additional selection criteria in declination, 5 < 5°, in order to observe galaxies from LCO. We 
then randomly sampled from the resulting catalog to generate a sample which would be observable 
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during the one year of allotted telescope time. Finally, we restricted the sample to two ranges in HI 
mass, selecting galaxies with either log[M#j /M ] < 9.0 (sample A) or \og[Mui /M & ] > 9.5 (sample 
B). This last selection criterion allows us to capture LSBGs in the maximum range of evolutionary 
states (i.e. high and low gas-to-star ratios) and total masses. 

There are 112 galaxies from the Impey et al. (1996) catalogue satisfying the above constraints. 
Of these, 65 galaxies were observed with the 40-inch telescope and 47 galaxies with the 100-inch 
telescope. 

In order to have an optimal wavelength baseline, we observed all galaxies through both the 
J and K s near- infrared filters (the latter defined by Persson et al. (1998)). Due to the critical 
importance of accurate surface brightness measurements, and given the low surface brightness 
of these galaxies, we chose not to use the H band. The strong and variable OH emission which 
dominates the sky background at H is not easy to remove, and could significantly affect the apparent 
surface brightness of the galaxies. In addition, prominent fringing appears in the H band, putting 
additional hurdles to deriving the correct surface brightness flux. The J band, on the other hand, 
is stable during time scales of ~ 20 minutes, and the K s bands is stable enough typically for time 
scales shorter than 5 minutes, when the outside temperature is below 15°C (Galaz 2000). Unlike 
many earlier near-IR studies of nearby galaxies, we have sufficient sky background in our images 
to do a reasonable job of sky subtraction. Most galaxies are much smaller than the ~ 2 x 2 arcmin 
FOV of both telescopes, with just a few having extents larger than 1.5 arcmin in the B band. 



3. Observations and Reductions 

All galaxies were observed by GG and EZ during photometric nights at LCO, using the 40- inch 
Swope telescope and the 100-inch du Pont telescope. Seven runs of 3 to 7 nights were allocated 
for this project. Observations started on February 23rd 1999 and finished on October 23rd 1999. 
A NICMOS3 HgCdTe array was used in both telescopes, but with different pixel scales: 0.60 
arcsec/pix at the 40-inch, and 0.348 arcsec/pix at the 100-inch. The corresponding FOV were 
2.56' x 2.56' and 1.48' x 1.48'. Observation procedures were the same as used in Galaz (2000), 
and consist of a sequence of short dithered exposures for all the scientific targets, including the 
standard stars. Due to the low surface brightness which characterizes these galaxies, exposure times 
were calculated based on the known blue surface brightness, extrapolated to the J and K s band, 
and corrected for the additional contribution of the thermal noise and sky emission. Typical total 
exposure times were between 20 min and 45 min in J and between 40 min and 75 min in K s . 

Considering the low surface brightness of galaxies, special care was taken with flat-fielding, in 
order to avoid undesirable artificial spatial fluctuations. A large series of dome-flats were taken for 
each run, as well as sky-flats every sunset. Additionally, the raw images of each night were used 
to construct super-flats, by combining dithered images, using a median sky value obtained after 
a sigma clipping algorithm to discard cosmic rays, bad pixels, and any residual contribution from 
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the galaxies themselves. Usually, no fewer than 80 images per filter were used to construct nightly 
super-flats, and only images where the galaxies were small compared the chip size were included 
in the image stack. Afterward, dome-flats, sky-flats and super-flats were compared. The largest 
differences between sky and dome-flats were ~ 1.5%. 

For image reduction, we used a modified version of DIMSUM 3 . The reduction procedure is 
completely analogous to that described in Galaz (2000), to which the reader can refer for details. 
Out of the 112 original selected LSBGs, we successfully observed 88; 39 from sample A (low HI 
mass) and 49 from sample B (high HI mass). 



In order to derive color information we calibrated our magnitudes system using the stan- 
dard star system defined by Persson et al. (1998). For the standard star photometry IRAF's 4 
"DAOPHOT" software was used. We adopted a 10 arcsec aperture radius for all the standard 
stars, which corresponds to the same aperture used by Persson et al. (1998). Typical seeing was 
1.0 arcsec. No less than 6 standards were measured each night, which allow us to derive calibrated 
J and K s magnitudes. The transformation equations were 



where capitals denote standard magnitudes, small characters denote instrumental magnitudes, a\ 
and b\ are the zero points, a 2 and b 2 are the color terms, and £j and ^ s are the extinction coefficients 
in J and in K s respectively. X is the mean airmass of the observation. 

We have computed the transformation parameters of equations (1) and (2) using "photcal", 
included in IRAF. When measuring the above coefficients, color terms and zero points were ini- 
tially computed independently for each night, verifying that they changed by only a small amount 
(less than 5%) from night to night; because these terms depend only upon the telescope/mirror- 
instrument-filter combination, they should be constant over each run. We then adopted a single 
color term and zero point for the entire run (based upon the night-to-night average), and re- 
computed the extinction coefficients for each night separately. The values obtained were very close 
to the expected values £j = 0.10 and ^ s = 0.08 in (Persson et al. 1998). The resulting average 
values typical of all our runs are shown in Table 1. The typical color term is relatively small, and 
given the observed values of J — K s (between 0.10 and 0.80), the contribution of this term is ~ 0.15 



3 DIMSUM is the Deep Infrared Mosaicing Software package developed by Peter Eisenhardt, Mark Dickinson, 
Adam Stanford, and John Ward. 

4 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



3.1. 



Calibrations 



J = 3 + ai -X£j + a 2 (j - k s ) 
K s = k s + b 1 -X( tks +b 2 {j-k s ) 



(1) 
(2) 
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in J, and therefore may be marginally important when computing the spatial distribution of colors 
through a given galaxy. 

3.2. Galaxy Photometry 

Given the heterogeneous luminosity profiles of the LSBGs, several different techniques were 
initially explored to optimize measurement of the apparent magnitudes: (1) circular apparent 
magnitudes, using DAOPHOT and the growth curve method, (2) total magnitudes within elliptical 
apertures, computed using SExtractor (Bertin and Arnouts 1996). SExtractor also computes 
isophotal areas which allow one to compute mean surface brightness within a specified isophote level. 
A comparison between these two techniques allows us to conclude that a more robust measurement 
of the flux is obtained with SExtractor magnitudes, given the elongated shape of many of the 
LSBGs in the sample (see Figures 2 and 3). In LSBGs with small ellipticities, DAOPHOT and 
SExtractor magnitudes differ by ~ 8%. In what follows, we use only magnitudes computed with 
SExtractor on elliptical apertures. Magnitudes were determined within an isophotal level set by 
the quality of the data, at a 3a level above the sky; the resulting sizes of the elliptical apertures are 
listed in Table 2. The major uncertainty arises from the sky level, with a mode of ~ 5% in J and 
~ 8% in K s . We have attempted to compare our derived magnitudes with those from the 2MASS 
survey, but found no overlap in the galaxies studied. 

Table 2 shows the J and K s magnitude and color, for the observed LSBGs. The heliocentric 
radial velocity from the 21-cm line, the size of the elliptical aperture used for photometry, and the 
HI mass content from Impey et al. (1996) are also shown. The absolute magnitudes M in Table 2 
were computed as 

- 25 (3) 

where m is the apparent magnitude, H is the Hubble constant expressed in km sec -1 Mpc -1 (we 
use Hq = 75 km sec -1 Mpc -1 in this paper), and v pec is the Virgocentric peculiar velocity in km/s. 
Throughout this paper we use v pec = 600 km/s, as quoted by Kraan-Korteweg (1986); Giovanelli et 
al. (1998). This corrective term is important for galaxies with cz < 2, 000 km/s. It is worth noting 
that this correction does not depend on passband, and therefore does not affect colors. Equation 
3 is a good approximation up to a recessional velocity cz ~ 25,000 km sec -1 , where cosmological 
corrections become important and the parameter go must be taken into account (from Table 2 we 
see that our largest velocity is ~ 23, 000 km sec -1 ). In equation 3 we have neglected /c-corrections. 
At the LSBGs redshifts /c-corrections are smaller than our photometric errors, and therefore are 
negligible. Note also that heliocentric radial velocity is measured from the HI emission (Table 2), 
slightly different from the optical heliocentric radial velocity. However, differences are smaller than 
100 km/sec, and therefore changes in absolute magnitudes are no larger than 0.05 mag. 

As can be seen from Table 2, LSBGs do not necessarily have faint absolute magnitudes in 
the near-IR. Instead, some of them appear to be as bright in the near-IR as normal high surface 



M = m — 5 los 



(cz 



J pec) 
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brightness galaxies. This immediately suggests that they do have a significant population of old 
stars. 



3.3. Measuring Central Surface Brightnesses and Scale Lengths 

Measurements of surface brightness and scale lengths give critical information on the light 
density and physical size of galaxies. We now use surface brightness fitting to measure the central 
surface brightness (as an indicator of the central light concentration) and scale length of both the 
disk and bulge components of the galaxies in our sample. The reader must recall that the surface 
brightness difference we measure for the same extended source located at different redshifts z\ < Z2 
is given by 

Af ,= 101og(I±f), (4) 

where Afi is expressed in mag arcsec -2 . In our case, our smallest radial velocity is ~850 km sec -1 
(z ~ 0.003), and our farthest galaxy is located at ~ 29,000 km sec" 1 (z ~ 0.1). This implies that 
the surface brightness difference for an object with a given surface brightness at z\, but measured at 
Z2 is, at most, A/j, = 0.4 mag arcsec -2 . However, at ±3<r level around the average redshift, we have 
v\ ~ 5,000 km sec -1 and V2 ~ 13,000 km sec -1 , which yields < A/j, >~ 0.08 mag arcsec -2 , which 
is inside the magnitude error bars. Therefore our surface magnitude estimates are not strongly 
modified by the range of distance included in the sample. 

We adopt the same approach as given by Beijersbergen, de Blok and van der Hulst (1999), 
and fit our data to the following parametric profile 

J(r) = / exp(-r//i), (5) 

or in terms of magnitudes per square arc-second 

fi(r) = + 1-086 x (r/h). (6) 

In equations 5 and 6 Iq and hq are the central surface brightness, the latter measured in mag 
arcsec -2 . The parameter h is the scale length, measured in arcsec and transformed to physical 
units (kpc). Similar to Beijersbergen, de Blok and van der Hulst (1999), we initially fit bulges 
using de Vaucouleurs profiles (de Vaucouleurs 1948). However, the profile was not a good match 
to the bulge light profile, and poor fitting results were obtained. Following Beijersbergen, de Blok 
and van der Hulst (1999), we then switched to fitting the bulge with an exponential light profile 
as well, yielding much better fits. An example bulge and disk fit is presented in Figure 1. The 
double exponential profile provides an adequate fit to the data over a large range in radius, though 
slight deviations suggest that an intermediate Sersic profile for the bulge might provide a slightly 
improved fit. We find that most of the galaxies require a bulge component in the near-IR (i.e. they 
show substantial curvature in the inner parts of a \i vs r plot), but some do not and represent 
completely disk dominated galaxies. 
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4. Results and Discussion 

4.1. Morphologies and Surface Brightnesses 

A mosaic with reduced and combined J images is shown in Figures 2 (sample A, low HI mass 
LSBGs) and 3 (sample B, high HI mass LSBGs). Typical image sizes are ~ 1.5' x 1.5' arcmin. Both 
Figures 2 and 3 and our surface brightness decompositions reveal a number of significant results 
for the morphologies and surface brightnesses of our sample. 

First, many of the LSBGs classified as late types based upon the optical B band image instead 
show a dramatic bulge in the near-IR (confirming the presence of prominent near-IR bulges noted 
by Beijersbergen, de Blok and van der Hulst (1999)). This result differs from that obtained by 
Bergvall et al. (1999) and Bell et al. (2000), who found similar morphologies in the B and J bands, 
but for a much smaller sample of galaxies. Examples of galaxies with different optical and near-IR 
morphologies (such as in the presence of a bulge in the near-IR) are the cases of LSBG #100 and 
LSBG #336 (for example, their J images in Figure 2). The fact that some LSBGs have high surface 
brightness bulges (fiQ(J) < 17.5 mag arcsec -2 ) has been commented by Beijersbergen, de Blok and 
van der Hulst (1999), who noted that many LSBGs could be regarded as high surface brightness 
bulges embedded in low surface brightness disks. Our results support such a hypothesis. In other 
words, while some galaxies have been classified as LSBGs on the basis of their disk properties, they 
are not low surface brightness overall. This difference is particularly striking in the near-IR, where 
the surface brightness more closely reflects the stellar surface density. Moreover, as discussed by 
Impey et al. (1996), many of the galaxies in the catalog from which our sample are drawn are not 
true LSBGs at all, and are instead indistinguishable from normal spirals in the optical. We find a 
similar result for our near-IR data. 

Second, in Figure 4 we plot the infrared luminosity of the galaxies vs the central surface 
brightness of the disk component. The figure shows that there is no correlation between these 
quantities and that almost all of the absolute magnitude vs surface brightness parameter space is 
covered. The upper left part of the distribution (large, luminous but LSB galaxies) is as populated 
as the the lower right part of the distribution (small, quite compact galaxies) . In our sample there 
is no bias against large and faint, nor small but luminous galaxies. We also observe in Figure 4 
a sharp cutoff towards bright surface brightnesses (at /j,j(disk) ~ 17.5 mag arcsec -2 ). This cutoff 
represents a real physical cutoff (i.e. nature does not make disks with higher surface brightnesses), 
rather than a selection effect propagating from an optical surface brightness cutoff in the B band. 
If it were the latter, then we would not expect the surface brightness of the cutoff to be constant 
with the near-IR absolute magnitude, given the strong dependence of the optical to near-IR color 
with galaxy luminosity. This physical cutoff further supports our belief that our sample does 
contain examples of galaxies with normal disk surface brightnesses, in addition to examples with 
disk surface brightnesses over 3.5 magnitudes fainter, thus giving us a true continuum of galaxies 
from normal disks to extremely low surface brightness systems. 
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Third, we have also included in Figure 4 a plot of the dependence of the bulge central surface 
brightness on the near-IR absolute magnitude. Unlike for the disks, we do see a correlation when 
overall surface brightness is considered (with a Spearman test suggesting that less than a 10% 
chance that they are uncorrelated for the low HI mass subsample). In addition, while Figure 4 
shows a cutoff at bright surface brightnesses for the disk component, the bulges in these same 
galaxies extend to much higher surface brightnesses (with a possible, but less sharp cutoff at 
H j (bulge) ~ 16 mag arcsec~ 2 ). For a comparison with the optical see Figure 4a from O'Neil, 
Bothun and Schombert (2000). 

Fourth, we compare the sizes (in kpc) and surface brightnesses of the bulge (left) and disk 
(right) components in Figure 5. It is apparent that LSBGs of the high HI mass subsample B (open 
circles) have brighter bulge surface brightnesses but similar bulge sizes, compared with LSBGs of 
the low HI mass subsample A (filled circles). A different trend, however, is observed for disks 
(right panel). In this case, galaxies of sample A have systematically smaller disk scale lengths 
with respect to those of sample B, but with little change in typical surface brightness. Therefore, 
galaxies of sample B have large disks and very high surface brightness bulges. Those galaxies with 
lower HI masses have small disks, and bulges that are structurally similar to those of the high HI 
mass sample, but that have very low surface brightness. The same behavior can be seen in the K s 
band (however we have better signal-to-noise in J). 

Rather than just comparing our two samples against each other, we can also compare mor- 
phology with HI content directly, by examining the bulge and disk structural parameters derived 
above as a function of HI mass. In Figure 6 the relationship is shown for the bulge and for the disk 
central surface brightness, in the left and the right hand panels, respectively. It is clear that the 
galaxies with low HI masses are almost entirely lacking in high surface brightness bulges, leading to 
a much smaller dispersion among the bulge central surface brightnesses. In contrast to the trends 
in bulge central surface brightness, there are no strong trends seen in the disk surface brightnesses. 
Instead, there is a very large dispersion among the disk properties, independently of the gas mass 
of the galaxies. 

Note that in all the above analyses using the near-IR bands, we are sampling nearly the same 
stellar populations both in the bulge and in the disk, namely red giants and old stellar populations. 
Thus, the relative sizes and surface brightnesses reflect true differences in the distribution of stellar 
density. Note also that these measurements of the current stellar content are not a priori related 
to the remaining HI content after star formation. However, the strong trends seen in Figure 5 and 
Figure 6 strongly suggest that the stellar surface density of a galaxy is indeed tightly coupled to 
the amount of neutral gas the galaxy currently has, though probably indirectly through the total 
galaxy mass, as we now demonstrate below. 
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4.2. Comparison between HI content and near-IR properties 

In the above analyses we have begun to consider how some properties vary as a function of 
the total HI mass. However, it is possible (if not likely) that the HI mass is not the sole controlling 
quantity driving these trends. It is quite possible that the HI mass (Mhi) correlates closely with 
the overall galaxy mass - i.e. that "bigger" galaxies have more of everything. 

In order to begin examining this possibility, Figure 7 shows the strong correlation we observe 
between the HI mass and the near-IR absolute magnitudes. Indeed, to first order, galaxies which 
have more stars (as measured by the near-IR luminosity) have more HI. Note that there are ~ 8 
magnitudes between the most HI massive galaxies and the less massive ones. Recall the striking 
fact that this correlation (which is well defined in both J and K s ) is between the mass of the HI 
gas, and the near-IR luminosity of the stellar populations; a priori there should not be a correlation 
between the mass in stars a galaxy has formed and the amount of gas it has left. 

In Figure 7, straight lines represent constant values of the ratio of gas mass to stellar luminosity 
for M#//L = 0.01 (long-dashed line), 0.1 (solid line), 1.0 (dotted line), and 10.0 (dashed line). 
Note that most LSBGs of our total sample have Mhi/L < 1.0 both in J and K s , with an average 
of Mhi/L ~ 0.1 and large scatter. While the scatter is over an order of magnitude in relative 
proportion of stars and gas, it is much smaller than the overall range of 8 orders of magnitude in 
overall gas mass, and thus does not wipe out the correlation between HI mass and total stellar 
luminosity. The overall trend in Figure 7 suggests that the division between our samples A and B 
is largely one mass, albeit with a wide range of gas mass to stellar luminosity ratios. 

We may further examine our sample for trends in how the gas mass to stellar luminosity ratio 
varies as a function of the overall galaxy luminosity. Figure 8 shows the logarithmic gas mass to 
stellar luminosity ratio as a function of J and K s absolute magnitudes (where we have adopted 
Mq(J) = 4.03 and M Q (K S ) = 3.33). Although there is considerable scatter, the lower luminosity 
galaxies have a higher fractional gas content, on average. The straight lines in Figure 8 are least 
square fits, giving slopes of 2.86 ± 0.46 in J and 3.35 ± 0.44 in K s . 

The gas mass to stellar luminosity ratio as a function of the absolute magnitude (Figure 8) 
indicates the relative importance of gas and stars in the total baryonic mass budget (albeit indi- 
rectly) . Galaxies with high Mhi /L in the near-IR have been historically less efficient at converting 
their gas into stars, due to either young mean age (i.e. a late start), low star-formation efficiency, 
large late-time gas infall, fewer past episodes of Hi-stripping, or some combination of the above. 
The fact that we see some correlation between Mhi/L and the galaxies' luminosities suggest that 
overall there are indeed trends between galaxies' star formation histories and their stellar mass. 
However, the large scatter suggests that it is not a perfect one-to-one relationship. We will ad- 
dress these issues in more detail in subsequent sections, where we attempt to reduce the scatter by 
calculating the actual stellar mass, rather than using the near-IR luminosity as a substitute. 

Another possible quantity which correlates with Mhi is the galaxy morphology. The study 
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of morphological types shows a modest correlation with HI content. Within our sample, a set of 
LSBGs with similar HI masses have morphological types spanning the entire Hubble sequence from 
dim to Sa. However, despite the large scatter, a mean trend does exist for morphological types 
over the whole HI mass interval. Galaxies of sample A (with low HI masses and, via Figure 7, low 
stellar luminosities) exhibit morphological types of dim, Sm, Sd, and only a few earlier type Sc 
and Sb. Only three cases are Sa. In contrast, LSBGs of sample B are much earlier type, and are 
classified in general as Sc, Sb and Sa, although some cases of Sm are also present (see Figures 2 
and 3, and text below). 

Much stronger trends are seen between the morphology of the galaxies and their gas mass to 
stellar luminosity ratio. Figures 9 and 10 show images of the galaxies sorted by their Mhi/Lk s 
values. The range in Mhi/Lk s goes from Mhi/Lk s = 0.028 - 0.128 M & /L & in Figure 9 (from 
LSBG #036 with \og(M HI /M Q ) = 8.74, to LSBG # 437 with log(M HI /M G ) = 7.92), and continues 
up to M HI /L Ks = 1.475 M Q /L Q in Figure 10 (galaxy LSBG #446 with log(M HI /M Q ) = 9.63). 
It is apparent in these two Figures that as Mhi/Lk s becomes larger, galaxies extend to much 
later Hubble types, beyond the text-book Hubble sequence. Galaxies with small Mhi/Lk s are 
earlier Sa/Sb/Sc types, while galaxies with large Mhi/Lk s tend to be Irr galaxies, suffer from tidal 
distortions, plumes and other peculiarities. 



4.3. Stellar mass-luminosity ratios 

In previous sections we have considered trends as a function of the HI mass to stellar luminosity. 
This ratio gives a measure of the evolutionary state of the galaxies, indicating how far the galaxy 
has proceeded along the path of converting gas into stars. However, the infrared stellar luminosity 
is only an approximate indicator of the stellar mass. We may obtain clearer and more physically 
meaningful results by using stellar population synthesis codes to translate the stellar luminosities 
into estimates of the actual stellar masses, allowing us to make a direct comparison between the 
gaseous and stellar mass. 

To convert from IR luminosity to stellar mass, we may take advantage of the work of Bell 
and de Jong (2001), which demonstrated that although the stellar mass-to-light ratio varies by up 
to a factor of 2 in the near-IR, there are clear trends between the stellar M/L and color. Using 
spectrophotometric spiral galaxy evolution codes, Bell and de Jong (2001) showed that one can 
use the colors of galaxies to derive their relative stellar mass-to-light ratios (under the assumption 
of a universal spiral galaxy initial mass function) . They also showed that uncertainties in the dust- 
reddening estimates do not strongly affect the final derived stellar masses of a stellar population. 

For testing whether the variations of M/L with J — K s color are robust to different stellar 
population synthesis models, we have used Bell and de Jong (2001) 's relations between the stellar 
M/L and the V — J and V — K colors, to derive the relationship between M/L and our J — K s 
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color 5 . Three different models are tested, which are (1) an instantaneous burst with a Scalo IMF, 
from recent Bruzual & Chariot models (Liu, Chariot and Graham 2000), (2) a recent version of 
PEGASE 6 models (Fioc and Rocca-Volmerange 1997, 2001) with an IMF slope of -1.85, and 
a model with a formation epoch with bursts, adopting a Salpeter IMF. All models have solar 
metallicity. The key point is that M/L vs color slopes do not depend on the IMF, which only 
affects the zero point. The stellar masses do not vary more than 7% between the two models. 
We have also used sub-solar metallicities and again the variation in the stellar masses are ~ 10%, 
and consistent with the variations between models. In what follows, we adopt the Salpeter model 
(using Table 1 from Bell and de Jong (2001)) and obtain coefficients for M/L as a function of 
J — K s color 7 . With the adopted relationship between stellar mass to light ratio and J — K s 
color, we can compute the actual baryonic masses Mb aryon i c from the derived stellar masses M ste u ar 
(= Lk s x (M/Lk)) and the HI mass as Mb ar yonic = M ste u ar + Mhi- We are implicitly assuming 
that the mass contribution from molecular, ionized, and metallic gas is negligible for our sample. 

Figure 11 shows the gas-to-star mass fraction as a function of stellar mass, after computing the 
stellar masses as described above. Figure 11 clearly shows that LSBGs with larger stellar masses 
tend to have smaller gas fractions, and LSBGs with small star masses have larger gas fractions. 

We have also included in Figure 11 diagonal lines representing constant HI mass; dashed line 
for Mhi = 7.0 M (close to our lowest HI mass), solid line for 9.0 M (upper limit of our sample 
A), dotted line for 9.5 M (lower limit for our sample B), and long-dashed line for 10.0 M (close to 
our largest HI mass). These lines can be used for assessing the importance of selection criteria. For 
example, there seems to be a lower limit for the mass fraction Hl-to-stars at ~ 5 x 10~ 2 . However, 
these LSBGs can exist but their HI would not be detectable, and thus they may not be present in 
the subsample of galaxies with measured HI fluxes that we selected from the Impey et al. (1996) 
catalogue. 

Figure 11 also confirms our conclusions (1) that low HI mass does not imply that a galaxy has 
"used up" its gas reservoir and converted it entirely to stars (given the existence of galaxies with 
low Mhi an d high Mhi /M s t ar ) and (2) that likewise high HI mass does not imply that a galaxy has 
only begun to convert HI into stars (given the large numbers of galaxies with large HI masses but 
large stellar masses as well). In other words, both high HI mass galaxies and low HI mass galaxies 
have a large range of stellar masses (with low HI mass galaxies having stellar masses between 10 7 
and 10 10 M , and high HI mass galaxies having stellar masses between 10 9 and 10 12 M ), and thus 
selecting on HI alone is not sufficient to identify galaxies which are gas-rich or gas-poor. 

We have also plotted similar relations as in Figure 11, but using the ratio of the HI mass to 

5 Differences between K and K s are less than ~ 0.08 mags. 
6 PEGASE is "Projet d'Etude des GAlaxies par Sinthese Evolutive. 

7 Bell and de Jong (2001) note that these relations have been calibrated for maximal disks, and may be slight 
overestimates if the disks are submaximal. 
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the total baryonic mass (M sta rs + Mhi), shown in Figure 12. Note that the ratio Mhi /Mb ary onic 
has an upper limit, when all the baryonic mass is concentrated in the gas (ratio = 1). We see 
again that increasing baryonic mass implies that a galaxy has converted a much larger fraction 
of its baryonic mass into stars. We have also included in Figure 12 the locus separating galaxies 
according to their dominant component, stars or gas. Again, note that "star-dominated" LSBGs 
can still have large quantities of HI, and similarly, "Hi-dominated" LSBGs can have a small Hi-mass 
if the overall baryonic mass is small. The extreme galaxy at Mbaryonic ~ 3 x 10 11 M Q is LSBG 224 
which, even with a large HI content (Mhi = 8-0 x 10 10 M Q ), has a baryonic mass represented almost 
exclusively by the stars. It is clearly a prominent Sa galaxy (see Figure 9). In general, all of the 
star-dominated LSBGs are Sa/Sb galaxies, whereas the Hi-dominated LSBGs have morphologies 
ranging from spirals to distorted or fuzzy irregulars. Among the most gas rich systems we find 
the most irregular and fuzzy LSBGs (Hi-rich galaxies) in our sample. Note that, although the 
number of galaxies from sample A and B is roughly similar, most of the galaxies lie below the 
dotted horizontal line, and therefore are star dominated ones. 

If we investigate these morphological trends further, we discover that both the stellar mass 
and the stellar-to-HI mass fraction are related closely to surface brightness. Figure 13 shows the 
bulge central surface brightness as a function of the ratio Mhi /Mbaryonic-, and as a function of the 
total baryonic mass. The ratio Mhi /Mb aryon i c increases with decreasing surface brightness, showing 
that at higher bulge surface brightnesses, galaxies have produced a much larger fraction of stars 
over their lifetime. It may be that the formation of a bulge creates an extremely efficient epoch 
of star formation, as gas is driven to very high central densities. Given the observed correlations 
between gas density and star formation rate (Kennicutt 1998), then we would expect that galaxies 
with the highest stellar densities (i.e. the high surface brightness bulges) had the highest past 
star formation rates, and thus the lowest current gas fractions. In other words, all galaxies are 
"born" at the locus of Mhi /Mbaryonic = 1> and then travel downwards and to the left in Figure 13 
to Mhi /Mbaryonic = 0, at a rate which depends in large part upon the central surface density of the 
gas. The end result is a lack of high surface brightness galaxies with large gas mass fractions, and 
low surface brightness galaxies with low gas mass fractions (e.g. the region around /j,q ~ 25 mag 
arcsec" 2 and M H i/M baryonic ~ 0.1). 

We may explore these ideas more directly by converting the observed surface brightnesses to 
stellar surface densities, using the color-dependent stellar mass to light ratios calculated above. In 
Figure 14 we have plotted the J — K s color and gas to baryonic mass ratio as a function of the 
central stellar surface density of the bulge S (in M /pc 2 ). We see a very strong trend (correlation 
coefficient of 0.22) towards lower gas fractions with increasing central surface densities, as expected 
for increased star formation efficiency at higher surface densities. The scatter is also reduced over 
Figure 13. In addition, we see a trend towards redder colors with increasing surface density, but 
we will defer a discussion of this effect until §4.4 below. 

While the increase in Mhi /Mbaryonic is a steady function of the bulge surface brightness, an 
equally striking trend is found with the the total baryonic mass. Figure 13 shows that the formation 
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of a high surface brightness bulge is a very sharp function of the overall baryonic mass. Galaxies with 
baryonic masses of less than Mb ar yonic = 10 10 M Q simply do not form bulges (using ^(bulge) < 18.0 
mag arcsec -2 , marked by a dash-dotted vertical line in Figure 13, as the demarcation of galaxies 
which do not have bulges significant enough to be detectable in the near-IR, i.e. those which cannot 
be distinguished clearly from the disk) 8 . 

The floor in the baryonic mass suggests that there is a tight coupling between the mass of the 
galaxy and the conditions necessary for bulge formation. We can envision such a correlation arising 
as follows. Assuming that there is no large scale segregation between baryons and dark matter, 
then the baryonic mass of the galaxy should be directly proportional to the total mass of the galaxy 
(M tota i /Mt, aryon ~ 16, for VL m = 0.25 and ^baryon = 0.015). Because the total mass of a galaxy is 
closely coupled to several astrophysical timescales for galaxy formation (such as the epoch of first 
collapse, or the end of an epoch of merging), one might expect to see a correlation between tracers 
of these quantities and the baryonic mass. In many scenarios, the bulge forms via a by-product of 
early merging in the galaxy, or via the collapse of low angular momentum material in the galaxy 
- both of which probably depend closely on the above timescales. Thus, it may be that future 
models of galaxy formation can use this observed floor for bulge formation as a critical test. 

The final question arises as to whether or not selection effects play a role in shaping the trends 
in Figure 13. For example, the lack of galaxies with high surface brightness bulges in the near-IR 
and low baryonic masses (i.e. the dashed oval) may be due to the fact that our sample was drawn 
from a catalog of low surface brightness galaxies. Such a catalog would be biased against bright 
bulges, unless accompanied by a low surface brightness disk, which would increase the overall mass 
of the galaxy, possibly pushing it up above the 1O 1O M floor. However, we suspect that the floor 
in Figure 13 is real, for several reasons. First of all, the existence of large numbers of galaxies 
with fiQ{bulge) > 18 mag arcsec -2 suggests that the sample by no means excludes galaxies of high 
near-IR surface brightness galaxies. Indeed, an analysis by Sprayberry et al. (1996) shows that 
the original Impey et al. (1996) sample includes a large fraction of high optical surface brightness 
galaxies, indistinguishable from "normal" spirals. Second, if the addition of a low surface brightness 
disk were necessary for inclusion in the original Impey catalog, there is no reason why the lower 
left of the plot should not be populated with low mass galaxies that would otherwise have dropped 
off the bottom of the plot. Finally, the observational selection effects would be most apparent in 
a plot of purely observational quantities (for example, Mj vs (io(bulge, J)), but in such a plot, we 
see the plane is evenly populated (e.g. Figure 4), and we can identify no steps in the process of 
converting from Mj to M-baryonic that would lead to the creation of an artificial floor in the baryonic 
mass associated with high surface brightness bulges. 

We stress that the fact that no gas poor but low surface brightness galaxies with very low 
values of Mm /Mb ar yonic and faint surface brightness appear in the Figure (/Kq ~ 25 mag arcsec -2 



8 For the galaxies with fj,o(bulge) > 18.0 mag arcsec 2 , we find mostly irregular morphologies, with a larger 
tendency toward showing spiral morphologies when dominated by stars (i.e. low Mijj/M stars ). 
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and y^Hl /^baryonic ~ 0.1), which could be a combined effect of our initial selection on surface 
brightness and gas content {hq(B) < 23.5 mag arcsec -2 ), even though we have not imposed a low 
cut for the HI mass. We know that large galaxies with bright total magnitudes but extremely low 
surface brightness and large disk scale lengths do exist, like Malin 1 (with Mb ~ —23.0, \xb = 26.5 
mag arcsec -2 and Mhi ~ 1O U M , Matthews, van Driel and Monnier-Ragaigne (2001)) and are 
placed at such a locus in the diagram, where spirals with faint bulges and large stellar masses lie. 
However, Figure 13 suggests that these galaxies are rare. 

4.4. The interplay between color and mass content 

In the previous section we have explored trends involving Mhi /Mb ar y 0nic as an indicator of the 
evolutionary state of a galaxy. If a galaxy steadily transforms its initial reservoir of gas into stars, 
then as M# j/M(, an/on j c decreases, the mean metallicity of the stellar population should increase as 
well. In this section we explore metallicity trends within our sample, after demonstrating that the 
near-IR J — K s color is an excellent indicator of the mean stellar metallicity. 

Within optical wavelengths, the variations in the broad band colors of galaxies are thought 
to be primarily driven by changes either in the mean age of the stellar population or in the dust 
content (Bell and de Jong 2001). However, in the near-IR the situation is quite different. First 
of all, reddening due to dust is negligible in all but the most extreme starburst galaxies. Secondly, 
and most importantly for this analysis, changes in the near-IR color of a galaxy are driven almost 
entirely by the mean metallicity of the galaxy; in the near-IR, the light from a galaxy is dominated 
by stars on the red giant branch, the location of which shifts to redder or bluer colors depending 
on metallicity. To demonstrate this effect, in Figure 15 we plot stellar population synthesis models 
from Bruzual & Chariot (Liu, Chariot and Graham 2000) of the J — K s color as a function of 
B — V, for different metallicities and exponentially declining or increasing star-formation rates (i.e. 
with different mean stellar ages), assuming that star-formation started 12 Gyr ago. The J — K s 
color index is clearly very metal sensitive, and relatively age insensitive; see Bell et al. (2000). At 
a single stellar metallicity, the change in J — K s with mean stellar age is never more than ~0.15, 
even when the mean age of the stellar population changes by 11 Gyr. This implies that whenever 
one witnesses a change of more than A(J — K s ) ~ 0.15 then one must be witnessing a significant 
change in stellar metallicity. We now apply these results to observations of the color distribution 
within our sample. 

In Figure 16 we show color-magnitude diagrams for galaxies appearing in Table 2. The plot 
clearly shows that LSBGs span a wide range in IR color (A(J — K s ) ~ 0.6), immediately suggesting 
a wide range in the stellar metallicities of the galaxies. As argued above, the largest possible shift 
in color due to age changes alone (A(J — K s ) ~ 0.15) is far too small to explain the full range of 
colors seen in our data. In other words, although changes in mean stellar age may contribute to 
the range in color, they cannot be responsible for all of it. 
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We may place some limits on the maximum range of metallicity compatible with the observed 
spread in J — K s color by assuming that the mean stellar age is constant (i.e. that all the color 
changes are due to metallicity alone). Based upon the grids in Figure 15, a factor of 10 change in 
the metallicity produces changes in color of roughly A( J — K S ) ~ 0.2 — 0.3. Therefore, our observed 
range in galaxy colors suggests a maximum change in mean stellar metallicity of more than a factor 
of 100 across the entire sample. Alternatively, we can calculate a minimum range in metallicity by 
allowing for the maximum spread in age (A(J — K s ) ~ 0.15). The remaining color difference still 
implies a minimum range of over a factor of 20 in mean stellar metallicity. These two limiting cases 
suggests that over the mass range in our sample, we have a range in stellar metallicity covering a 
factor of 20 - 100. 

In addition to demonstrating the overall range of near-IR colors, Figure 16 shows a statistically 
significant trend between the IR color and the stellar luminosity of the galaxy, suggesting a trend 
between metallicity and stellar mass. (We may also use Figure 16 to compare the typical colors of 
the low gas-mass and high gas-mass subsamples that the lower HI mass galaxies in sample A are 
bluer that sample B by ~ 0.3 magnitudes, on average, immediately suggesting a mass dependent 
trend in galaxy metallicity, given the trends between HI mass, total stellar mass, and baryonic mass 
between the two samples). We can explore these trends more explicitly in Figure 17, where we plot 
J — K s directly as a function of HI mass. As expected there is a systematic trend between the HI 
mass of a galaxy and its near-IR color. Finally we plot J — K s as a function of the stellar mass and 
of the total baryonic mass in Figure 18. Again, we see trends in the near-IR color with the stellar 
and baryonic masses (with a slope of 0.22 ± 0.03 and a rms of 0.12 for the linear fit of color versus 
logarithmic stellar mass, and with a slope of 0.27 ± 0.03 with a rms of 0.13 for the linear fit with 
logarithmic total baryonic mass). 

The general picture painted by Figures 16-18 is one where as the stellar mass, gas mass, and/or 
total baryonic mass within a sample of galaxies increases, the galaxies becomes redder in J — K s , 
and thus more metal rich on average. However, around all of these relations there is substantial 
scatter, suggesting that none of these quantities is the principal driver of metallicity. 

Both our intuition and our data suggest that the true driver of the metallicity is the gas mass 
fraction M# / /Mf, aryonic , which indicates the degree to which a galaxy has completed its lifetime of 
star formation. Galaxies that have only just begun to create stars have only begun to enrich their 
gas reservoirs, and thus should have a low mean stellar metallicity and very blue J — K s colors. 
Likewise, galaxies that have completely converted their gas into stars via several generations of 
star formation will have greatly enriched the gas from which the last half of the stars formed, 
leading to much more metal rich stellar populations and redder J — K s colors. In Figure 19 we 
plot the gas mass fraction Mhi /Mb ar yonic as a function of the metal sensitive J — K s color, and 
overall the behavior is just as expected. Compared to the previous figures, there is a substantial 
tightening of the relationship, especially towards the more stellar dominated systems (i.e. low values 
of Mni/Mbaryonic)- This suggests that the typical stellar metallicity of a galaxy is indeed largely 
set by the degree to which it has already converted its gas into stars. We have already seen hints of 
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this effect in Figure 14, where the J — K s color is highly correlated with the central surface density 
of the galaxy, which in turn sets the efficiency of converting gas into stars and thus the subsequent 
value of M HI /M baryonic . 

We may also gain additional insight by considering the scatter in Figure 19. In a "closed box" 
model for chemical enrichment, the mean metallicity of the stellar population increases monoton- 
ically with decreasing gas mass fraction, leading to a perfect one-to-one relation between the two 
quantities. If the closed box model is valid on the scale of galaxies, then the relationship between 
J — K s and Mhi lM baryonic should have zero scatter (modulo the small differences in J — K s pro- 
duced by changes in mean age). For the galaxies in Figure 19 where star formation has proceeded 
largely to completion (Mhi /Mb aryonic < 0.2), the scatter is comparable to our observational errors, 
suggesting that evolution could have proceeded roughly along a closed box pathway 9 , as also con- 
cluded by Bell and de Jong (2001) for normal spirals. In contrast, for more gas rich galaxies, the 
scatter in J — K s is significantly larger than our uncertainties, suggesting that the closed box model 
is not a good approximation to the gas and metallicity evolution of these systems. Instead, these 
systems may still be in a phase of sporadic gas infall, as opposed to the more massive star-dominated 
galaxies where the gas infall seems to have run to completion. Note also that these gas rich galax- 
ies span a much wider range in stellar luminosity (and thus stellar mass) at fixed Mhi /M barycmic , 
possibly suggesting a much wider range of mean stellar age, which could contribute to the larger 
scatter as well. 

To facilitate comparison with earlier work in the literature, we have also included in Figure 20 
plots of the J — K s color as a function of M^//L, which has traditionally been used as an indicator 
of the relative fraction of gas and stars. The observed trend toward lower Mhi/L and redder colors 
is observed in both bandpasses, and while it is more obvious in the K s band, the relationship is still 
much weaker overall than the trend with Mhi lM baryonic in Figure 19. We fit the trend with straight 
lines, giving a slope of —0.05 ± 0.04 in J (which is not significant), but a slope of —0.13 ± 0.04 in 
K s , which is. Note that O'Neil, Bothun and Schombert (2000) found no relationship between the 
gas fraction and the B — V color, due to the large scatter (in that work the gas fraction and color 
relationship arise from a compilation of different sources of data). Our observed tightening of the 
relationship reflects Bothun (1982) 's observations that near-IR colors are less scattered than blue 
colors when plotted against the gas fraction. 

We may explore the possible effects of age by considering previous work on the optical colors of 
LSBG's. For example, one may compare Figure 17 with Figure 5a presented by O'Neil, Bothun and 
Schombert (2000), where no correlation exists between the HI mass and the B — V (optical) color. 
The lack of any correlation between HI mass and the age sensitive B — V color index suggests 
a stochastic star formation history, as also argued by Schombert, McGaugh & Eder (2001) for 
LSB dwarfs. The measured colors are luminosity weighted, and thus very small changes in the 



The gas poor galaxies are also those with the most substantial bulges, and thus it may be that the closed box 
approximation is valid for the bulges alone. 
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current star formation rate produce large changes in the luminosity-weighted mean age, and thus 
in B — V . As quoted by O'Neil, Bothun and Schombert (2000), this clearly suggests the existence 
of "dormant" populations of LSBGs galaxies, which can have large quantities of HI but low current 
star formation. 

If we take the increased scatter in Figure 19 towards gas rich galaxies as an indication of 
continuing gas infall, then we may simultaneously have an explanation for the sporadic past star 
formation history suggested by the optical data, namely gas density dependent star formation 
thresholds. In this scenario, ongoing gas infall occasionally raises the gas surface density above a 
critical threshold, triggering a burst of star formation and a strong bluing of the optical colors, 
whiled simultaneously diluting the metallicity of the gas phase and of subsequent generations of 
stars. This would yield a larger scatter both in metallicity and in age, leading to the much larger 
spread in J — K s colors. It would also help to explain why many of the optically faint LSBGs in 
our sample are sufficiently old to have developed significant red giant branch populations, in spite 
of what appear to be young luminosity weighted ages implied by their typically blue optical colors. 

Finally, we note that all of the results above are based upon the color of the entire galaxy. If 
we were to consider only the disk component, the relationship between gas-to-star mass-luminosity 
ratio and color would be loosened with the increased scatter being due to the still low S/N which 
is obtained for the disk photometry in the near-IR. 

5. Conclusions 

The near-IR photometry of low surface brightness galaxies presented in this paper constitutes 
one of the largest database published until now. Because of our original selection of galaxies from 
the Impey et al. (1996) catalog, our sample includes not just a very large number of low surface 
brightness galaxies, but a number of normal spiral galaxies as well. 

The total sample studied provides good indication that an old, red population of low surface 
brightness galaxies does exist. Many of the LSBGs analyzed exhibit clear bulges and have high 
surface brightness disks in the near-IR, even though they were originally identified in a catalog of 
optical LSBGs 10 . Morphologically, the sample spans the full range of spiral Hubble types, from Sa 
to Im. While there are broad correlations between the galaxies' optical and near-IR morphologies, 
we find several cases of galaxies which appear to be late-type in the optical, but have prominent 
bulges in the near-IR (e.g. LSBGs # 100, 338, 384, and 446 in Figures 2 and 3). 

In addition to the morphologies, a quantitative analysis of the surface brightness profiles of 
the galaxies suggests that our sample contains galaxies with disk surface densities comparable to 
normal spirals, as well as a very high fraction of truly low surface brightness (and/or low surface 



10 These bulge dominated LSBGs also tend to be red, as found by Bell et al. (2000) for the optically selected red 
LSBGs of O'Neil, et al. (1997) 
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density) disks. On the bright end, we find that the central surface brightnesses of the disks in 
our sample show a well defined cutoff brighter than /j,j(disk) ~ 17.5 mag arcsec -2 (independent of 
the J-band absolute magnitude). We argue that this cutoff is physical, and reflects the maximum 
surface density at which disks are stable. On the LSBG end, our sample includes a very high 
fraction of galaxies that are indeed low surface brightness even in the near-IR, with some reaching 
roughly 3.5 magnitudes fainter in surface brightness than the bright cutoff. The disks also have 
sizes that are systematically smaller for lower mass galaxies. 

As part of our analysis, we have calculated the total baryonic masses (M;, an/on j c = M gas + 
M stars) of galaxies in our sample. For M gas , we have taken the HI mass from the Impey et 
al. (1996) catalog, and assumed that other gaseous components are negligible. For M sta rs, we have 
used the color-dependent stellar mass-to-light ratios of Bell and de Jong (2001) to translate the 
observed near-IR luminosity into stellar mass. The mass-to- light ratios also allow us to convert 
observed surface brightnesses into stellar mass surface densities. The resulting masses and surface 
densities allow us to consider trends involving fundamental physical quantities (such as the ratio 
of the gas mass to the total baryonic mass). The total baryonic mass itself is probably a good 
indicator of the total mass of the galaxies, assuming that baryon and dark matter are evenly mixed 
on large scales. 

Among the masses of the various components, we see strong correlations, many of which follow 
trends previously identified in the optical (for example, see compilations in Schombert, McGaugh &: 
Eder (2001)). Galaxies with large stellar masses tend to also have high HI masses, higher central 
surface densities, but small gas mass fractions (i.e. are Hi-poor). Likewise, galaxies with small 
stellar masses tend to have low HI masses but have a larger fraction of their baryonic mass in the 
gas phase (are Hi-rich, as for the LSBGs studied by McGaugh and de Blok (1997)). We note that 
even galaxies which have formed large masses of stars can still retain large reservoirs of neutral gas, 
even after rich episodes of star formation. 

We find several strong morphological correlations with the masses of the galaxies. There are 
systematic trends towards more irregular, diffuse morphologies and lower surface densities with 
increasing gas mass fraction (or alternatively, M^//L^ 3 ). We also find that the central surface 
brightness of the bulge component varies strongly with the near-IR luminosity and HI mass, as well 
as with the stellar and the total baryonic mass (such that the bulges have higher surface brightness 
for larger values of the listed quantities). Moreover, we find that our sample exhibits a "floor" in 
the baryonic mass required for bulge formation, at 1O 1O M ; we have no galaxies with lower baryonic 
masses that have formed high surface brightness bulges. 

We have also explored how the near-IR color varies within our sample. In order to interpret 
the colors, we have used stellar population synthesis models to argue that the J — K s color is highly 
metal sensitive, such that redder colors imply higher metallicities. We argue that based upon the 
range of observed colors our sample spans a factor of more than 20 in metallicity (but less than a 
factor of 100). We find that galaxies have redder J — K s colors with increasing HI mass, stellar 
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mass, and baryonic mass, suggesting that more massive galaxies have in general have more metal 
enriched stellar populations. 

We find very tightest correlations between color and the gas mass fraction Mhj /Mb ar yonic in 
particular. These correlations are significantly tighter than previous optically determined trends, 
even in the redder V — I color (see compilation in Figure 6 of Schombert, McGaugh & Eder (2001)). 
In a closed box model for chemical evolution, there should be a one-to-one relationship between 
the metallicity and the gas mass fraction. At small values of Mhi /^-baryonic: t ms is essentially 
what we see, in that the observed scatter in the metal-sensitive J — K s color is comparable to our 
observational uncertainties. This suggests that in the gas poor, bulge dominated massive galaxies, 
chemical evolution has proceeded along the canonical closed box pathway (possibly confined to the 
bulge, which dominates the luminosity and thus the colors of the galaxy). 

We find a different behavior among the gas rich, disk dominated, lower mass galaxies. In 
these systems there is large scatter around the relationship between color and Mhi /^baryonic- This 
suggests that not only is the enrichment of these galaxies on-going (given their high gas mass 
fraction), but that it is not proceding as expected for a "closed box" model of enrichment. Instead, 
it is likely that these gas rich systems are continuing to experience episodes of gas infall. These 
episodes of gas accretion first increase the ratio of gas to stars, and then dilute the metallicity of 
subsequent generations of stars. 

Episodic accretion would also explain the stochastic star formation suggested by analyses 
of the optical colors of LSBGs (O'Neil, Bothun and Schombert 2000; Bell and de Jong 2001) 
and theoretical calculations (Gerritsen and de Blok 1998; Jimenez et al. 1998). Temporary 
increases above a gas surface density threshold for star formation could dramatically increase the 
star formation rate (Kennicutt 1998), triggering bursts of star formation and an increased rate of 
gas consumption. This change in the luminosity-weighted mean stellar age would produce additional 
scatter in the color, scatter which will be more dramatic in the optical. 

Within our sample, we see strong evidence that the star formation efficiency is indeed surface 
density dependent. We find strong correlations between the gas mass fraction Mhi /M^ arycmic and 
the central stellar surface density (i.e. M stars /pc 2 ), suggesting that galaxies with the highest initial 
gas surface density have also had the highest efficiency in turning gas into stars. This link between 
surface density and star formation efficiency also suggests that the young mean stellar age of LSBGs 
(as deduced from luminosity- weighted optical colors) results principally from the galaxies being slow 
to convert gas into stars, rather than from their being slow to assemble (see also Bell et al. (2000), 
Schombert, McGaugh & Eder (2001), de Blok & van der Hulst (1998)). However, the evidence for 
on-going gas accretion in lower luminosity galaxies suggests that this second effect does play some 
role as well. 

For future work, we intend to extend the current work with the addition of optical colors from 
our on-going database of B and R photometry for the same galaxies presented here, as well as from 
other bands provided by The Sloan Digital Sky Survey. We shall compute dynamical masses, and 



-23- 



current gas phase metallicities using also future spectroscopic data. This will allow us to perform 
more critical and accurate tests of the ideas presented in this work. 

We acknowledge Las Campanas Observatory for granting the large investment in observing 
time necessary for this project. We also thank all of the staff at Las Campanas. We would like 
to acknowledge Pablo Araya in the reduction of some images appearing in this paper, as well as 
Andreas Reisenegger for comments and shaking discussions around the subject presented here. We 
also thank the anonymous referee for many helpful suggestions that greatly improved the paper. 
Part of this research was funded by the Andes-Carnegie fellowship C-12927, during a post-doctoral 
position at Las Campanas Observatory (GG). GG thanks funding from "Proyecto DIPUC 2001/14- 
E". GG and LI thanks "Proyecto FONDAP" Center for Astrophysics. LI acknowledges "Proyecto 
Puente PUC". JD was partially supported by NSF grant AST-990862 and by the Alfred P. Sloan 
Foundation during this work. 

This research has made use of the NASA/IPAC Extragalactic Database, which is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National 
Aeronautics and Space Administration. The research has also made use of the VizieR service from 
the Centre de Donnees Astronomiques de Strasbourg (CDS), and the Digitized Sky Survey, which 
is produced at the Space Telescope Science Institute under U.S. Government grant NAG W-2166. 
The images of these surveys are based on photographic data obtained using the Oschin Schmidt 
Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates were processed into 
the present compressed digital form with the permission of these institutions. 

REFERENCES 

Bergvall, N., Ronnback, J., Masegosa, J., Ostlin, G. 1999, A&A, 341, 697 
Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 

Beijersbergen, M., de Blok, W., k van der Hulst, J. 1999, A&A, 351, 903 
Bell, E., & de Jong, R. ApJ, 550, 212 

Bell, E., Barnaby, D., Bower, R. et al. 2000, MNRAS, 312, 470 
Bell, E. & de Jong, R., 2000, MNRAS, 312, 497 

Bell, E., Bower, R., de Jong, R., Hereld, M., & Rauscher, B. 1999, MNRAS, 302, 55 
de Blok, W., & McGaugh, S. 1997 MNRAS, 290, 533 
de Blok, W., & van der Hulst, J., A&A, 335, 421 

de Blok, W., McGaugh, S., & van der Hulst, J. 1996, MNRAS, 283, 18 



- 24 - 



de Blok, W., van der Hulst, J., k Bothun, D. 1995, MNRAS, 274, 235 
Bothun, G. 1982, ApJS, 50, 39 

Bothun, G., Impey, C, Malin, D., k Mould, J. 1987, AJ, 94, 23 
Bothun, G., Impey C, k McGaugh, S. 1997, PASP, 109, 745 
Courteau, S. 1996, ApJS, 103, 363 

Dalcanton, J., Spergel, D., k Summers, F. 1997, ApJ, 482, 659 
Disney, M. 1976, Nature, 263, 573 

Fioc, M. k Rocca-Volmerange, B. 1997, A&A, 326, 950 
Fioc, M., k Rocca-Volmerange, B., in preparation 
Fisher, J., k Tully, R. 1975, A&A, 44, 151 

Freeman, K. 1970, ApJ, 160, 811 (erratum 1970, ApJ, 161, 802) 

Galaz, G. 2000, AJ, 119, 2118 

Gerritsen, P. k de Blok, W. 1998, A&A, 342, 655 

Giovanelli, R., Haynes, M., Salzer, J., Wegner, G., da Costa, L., k Freudling, W. 1998, AJ, 116, 
2632 

Impey, C. k Bothun, G. 1989, ApJ, 341, 89 

Impey, C, Sprayberry D., Irwin M., k Bothun G. 1996, ApJS, 105, 209 
de Jong, R. 1996, A&A, 313, 45 

Jimenez, R., Padoan, P., Matteucci, F., Heavens, A. 1998, MNRAS, 299, 123 

Kennicutt, R. 1998, ARA&A 

Knezek, P., k Wroten, J. 1994, A AS, 185, 77.13 

Kraan-Korteweg, R. 1986, A&AS, 66, 255 

Liu, M., Chariot, S., k Graham, J. 2000, ApJ, 543, 644 

Longmore, A., Hawarden, T., Goss, W., Mebold, U., k Webster, B. 1982, MNRAS, 200, 325 
Matthews, L., van Driel, W., k Monnier-Ragaigne, D. 2001, A&A, 365, 1 
McGaugh, S., k de Blok, W. 1997, ApJ, 481, 689 



-25- 



McGaugh, S. 1994, ApJ, 426, 135 

McGaugh, S. 1992, PhD Thesis, Univ. Michigan 

O'Neil, K., Bothun, G., k Schombert, J. 2000, AJ, 119, 136 

O'Neil, K., Bothun, G., Schombert, J., Cornell, M., k Impey, C 1997, AJ, 114, 2448 
Persson, E., Murphy, D., Krzeminski, W., Roth, M., k Rieke, M. 1998, ApJ, 116, 2475 
Schelegel, D., Finkbeiner, P. k Davis, M. 1998, ApJ, 500, 525 
Schombert, J., McGaugh, S. S., k Eder, J. A. 2001, ApJ, 95, 1389 
Schombert, J. k Bothun, G. 1988, AJ, 95, 1389 

Sprayberry, D., Bernstein, G., Impey, C, k Bothun, G. 1995, ApJ, 438, 72 
Sprayberry, D., Impey, C., k Irwin, M. 1996, ApJ, 463, 535 

Tully, R., Pierce, M., Huang, J.-S., Saunders, W., Verheijen, M., k Witchalls, P. 1998, AJ, 115, 
2264 

Tully, B. k Verheijen, M. 1997, ApJ, 484, 145 

Tully, R. k Fisher, J. 1977, A&A, 54, 661 

van den Bergh, S. 1959, Pub. David Dunlap Obs., 2, 159 

van den Bergh, S. 1998, Galaxy Morphology and Classification (Cambridge University Press) 

van den Hoek, L., de Blok, W., van der Hulst, J., k de Jong, T. 2000, A&A, 357, 397 

van der Hulst, Skillman, E., Smith, T., Bothun, G., McGaugh, S., k de Blok, W. 1993, AJ, 106, 
548 

de Vaucouleurs, G. 1948, Ann. d'Astrophys. 11, 247 

Verheijen, M. 1997, Ph.D. thesis, Univ. Groningen 

Zwicky, F. 1957, Morphological Astronomy (Berlin: Springer) 



This preprint was prepared with the AAS IATprjX macros v5.0. 



-26- 



Figure captions 

Fig. 1. — Light profiles in the J and K s bands for one LSBGs studied in the sample B. The total 
light profile is fitted by two exponential profiles, a bulge and a disk components, using equation 6. 

Fig. 2. — Mosaic of sample A (low HI mass sample). The label is correlative with column 1 of 
Table 2. Figure under request. 

Fig. 3. — Mosaic of sample B (high HI mass sample). The label is correlative with column 1 of 
Table 2. Figure under request. 

Fig. 4. — J absolute magnitude as a function of the disk central surface brightness. Filled circles 
indicate LSBGs of sample A, and open circles LSBGs of sample B. 

Fig. 5. — Bulge (left) and disk (right) central surface magnitudes as function of the respective scale 
lengths in kiloparsecs. Filled circles indicate LSBGs with \og\Mni /Mq] < 9.0 (sample A), and 
open circles \og[Mni /M@] > 9.5 (sample B). 

Fig. 6. — Bulge (left) and disk (right) central surface brightness as a function of the HI mass. Note 
the large dispersion for a given HI mass. 

Fig. 7. — Absolute magnitude as a function of HI mass. Left panel is for J and right panel for K s . 
Diagonal lines represent constant gas fractions. Log(M#//L) = 0.01 for long-dashed line, 0.1 for 
solid line, 1.0 for dotted line, and 10.0 for dashed line. 

Fig. 8. — J and K s absolute magnitudes as a function of the logarithmic gas fraction. The straight 
lines represent least square fits, giving a slope of 2.86 ± 0.46 for J and 3.35 ± 0.44 for K s . 

Fig. 9. — J images of low surface brightness galaxies included in this paper, sorted by the gas 
fraction M HI /L K s. Galaxy 036 has a value M HI /h Ks = 0.028, and galaxy 437 M HI /L Ks = 0.128. 
Most of these LSBGs fall rather well in the Hubble sequence and most of them present clearly a 
bulge in near-IR passbands. Figure under request. 

Fig. 10. — J images of low surface brightness galaxies included in this paper, sorted by the gas 
fraction M HI /L Kg . Galaxy 225 has a value M HI /L Ks = 0.129, and galaxy 446 M HI /L Kg = 1.475. 
Most of these LSBGs are Irr types, but some of them are difficult to classify following the Hubble 
sequence. Most of these LSBGs do not present a bulge in the near-IR passbands. Figure under 
request. 

Fig. 11. — Hl-to-stellar mass ratio and color J — K s as a function of the computed stellar mass, 
which has been computed using results of Bell and de Jong (2001). See text for details and 
discussion. Diagonal lines represent constant gas mass dashed line for Mhi = 7.0 M ( close to 
our lowest HI mass), solid line for 9.0 M Q (upper limit of our sample A), dotted line for 9.5 M Q 
(lower limit for our sample B), and long-dashed line for 10.0 M Q (close to our largest HI mass). 
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Fig. 12. — The ratio Mhi /Mb aryonic as a function of the baryonic mass (HI mass + stellar mass). 
Diagonal lines represent constant HI content, as indicated in Figure 11. Dotted boxes indicate the 
locus of gas and star dominated LSBGs. 

Fig. 13. — Gas-to-stars mass fraction as a function of the bulge surface brightness. The labels 
indicate the locus where gas and star dominates and the their bulge brightness. The cutoff at 
Ho(J, bulge) = 18.0 mag arcsec -2 represent the surface magnitude limit where galaxies become to 
be either gas or star dominated. See text for discussion. 

Fig. 14. — The J — K s color and gas to baryonic mass ratio as a function of the bulge stellar 
surface density (measured in solar masses per sqr-parsec). The dashed lines represent linear fittings: 
log(M Hl /M baryonic ) = (0.40 ± 0.03) - (0.39 ± 0.06)1! and J — K s = (0.25 ± 0.02) + (0.19 ± 0.04)E. 
The correlation coefficient is 0.22 and 0.15 in the upper and lower panel, respectively. 

Fig. 15. — Color grid for the J — K s index as a function of the B — V index, for different stellar 
formation rates and metallicities. Metallicities range from 0.005Z Q to 2.5Z Q . Top labels denote 
different exponential star-formation rates, where oo denotes a constant star-formation rate. Bottom 
labels denote mean ages in Gyr. The star-formation started 12 Gyr ago for all bursts. 

Fig. 16. — Near-IR J (left) and K s (right) color-magnitude diagrams for LSBGs presented in 
this paper. Filled circles indicate galaxies with \og[Mni /Mq] < 9.0 (sample A), and open circles 
galaxies with log [Mr- j/M Q ] > 9.5 (sample B). 

Fig. 17. — J — K s colors as a function of the HI mass. There is a clear positive trend in the sense 
of larger Mhi having redder J — K s . This trend is not observed in optical colors, as discussed in 
the text. 

Fig. 18. — Computed stellar and baryonic masses as a function of the color J — K s . Straight lines 
represent linear fittings. Both slopes are very similar, being 0.23 ±0.03 for the stellar masses (solid 
line), and 0.27 ± 0.03 for baryonic masses (dashed line). 

Fig. 19. — Gas-to-baryonic mass fraction as a function of the J — K s color and K s absolute 
magnitude. Note that color become redder for smaller values of Mhi /Mb aryon i c . Dotted lines 
enclose the locus of star and HI dominated LSBGs. The scatter could reflect a star formation 
process using fresh HI during the evolution to redder colors, and so larger metallicity. See text for 
discussion. 

Fig. 20. — J — K s color as a function of the gas fraction computed in the J and K s bands. 
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Table 1. Average Zero points, color terms and extinction coefficients. 
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f 1 ) Average zero point for the J filter. 
(2) A 

vcragc color term for the J filter. 

Average zero point for the K s filter. 
^Average color term for the K s filter. 
( 5 ) Average extinction value for the J filter. 

Average extinction value for the K s filter. 

Number of nights considered for computing the average values. 
(8) Number of standard stars observed per night. 

standard deviation for the corresponding averages arc included between parenthesis. 
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Table 2. The sample of low surface brightness galaxies observed in the near-IR. 
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a Number correlative to the Impcy ct al. (1996) catalogue. 
k Name from the Impey et al. (1996) catalogue. 

c Heliocentric radial velocity in km/sec from the 21-cm HI line (Impey et al. 1996). 
d HI mass expressed as log(M H j/Mq ) (Impcy ct al. 1996). 
e Total exposure time in the corresponding filter. 

f Average between the major and minor axes of the elliptical aperture used 
by SExtractor, in arcsccs. The value corresponds to 3ct above sky, in 
the J filter. The same elliptical aperture is used for the K s image. 

^Calibrated aperture magnitude. 
^Photometric error. 

Absolute magnitude computed following equation 3, but considering v pec — 0. 
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